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iVe model the response of a state of the art micro-hole single-stage charge amplification device (‘microbulk’ Micromegas) in a 
gaseous atmosphere consisting of Xenon/trimethylamine at various concentrations and pressures. The amplifying structure, made 
vith photo-lithographic techniques similar to those followed in the fabrication of gas electron multipliers (GEMs), consisted of a 
LOO /rm-side equilateral-triangle pattern with 50 /rm-diameter holes placed at its vertexes. Once the primary electrons are guided 
nto the holes by virtue of an optimized field configuration, avalanches develop along the 50 ^m-height channels etched out of the 
rriginal doubly copper-clad polyimide foil. In order to properly account for the strong field gradients at the holes’ entrance as well 
is for the fluctuations of the avalanche process (that ultimately determine the achievable energy resolution), we abandoned the 
lydrodynamic framework, resorting to a purely microscopic description of the electron trajectories as obtained from elementary 
:ross-sections. We show that achieving a satisfactory description needs additional assumptions about atom-molecule (Penning) 
transfer reactions and charge recombination to be made. 
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. Introduction 

Microbulk Micromegas, introduced in [1], are a new 
eneration of charge-amplifying Micromegas devices [2] 
capable of delivering, on a selected mm^ spot, intrin¬ 
sic energy resolutions down to 10.5%@5.9keV (Neon/i- 
C4H10), 11.2%@5.9keV (Argon/i-C4Hio) and 7.3%@22keV 
(Xenon/TMA) at atmospheric pressure [1,3,4]. A similar 
hierarchy (He/Ne/Ar/Xe) has been experimentally ob¬ 
served for micro-meshes under sub-100 /rm^ irradiation in 
[5], and it can be linked to the efficiency of the avalanche 
multiplication process, as theoretically anticipated [6]. 

It is indeed remarkable that the performance of microb¬ 
ulk readouts deteriorates only slightly on real-sized exper¬ 
imental systems, showing 8.5%@30keV (Xenon/TMA) on 
700 cm^ [7], and 14.7%@5.9keV (Argon/i-C4Hio) on 15 
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cm^ [8]. Furthermore, the 730 normal liter TPC in [7] was 
operated steadily under irradiation from a 15 kBq 7-ray 
source for more than 100 live days. A modest x2 reduction 
of the maximum working gain, relative to values obtained 
earlier for small 3.5 cm^ wafers in [4] was observed. Given 
the large number of holes typically employed (^ 10®/m^), 
such good scaling properties require a great accuracy on 
the hole diameter and gas gap, at the %-level or below. 

Besides offering an accurate mechanical construction, ca¬ 
pable of delivering a good energy resolution and stability 
on large areas, additional benefits of the microbulk tech¬ 
nique are the low material budget and relatively high ra¬ 
diopurity [9]. Their overall performance makes them attrac¬ 
tive for the construction of large time projection chambers 
(TPCs) for experiments on ‘Rare Event Searches’ under a 
broad range of admixtures and pressures. For some appli¬ 
cations, chiefly solar axion detection, radioactivity levels of 
10“^/s/cm^/keV [8] could be achieved in the sub-10 keV 
energy region for underground operation. 





Microbulk Micromegas are presently used in CAST [8] 
and n-TOF [10] experiments, and its application is foreseen 
or considered in a new generation of experiments involv¬ 
ing axions (lAXO [11]), dark matter (TREX-DM [12]) and 
PPOv decay [13]. The latter, naturally focusing on high- 
pressure ^^®Xe, builds on the pioneering works on charge- 
readouts performed by the Gotthard collaboration [14] and 
more recently in [7], and will try to challenge the state of 
the art of PPOv detection in gas, based on light-readout and 
presently led by the NEXT experiment [15,16]. 

Outside the ideal conditions experimentally realized in 
[5], for the practical operation of Micromegas readouts 
an interplay exists between geometrical tolerance, opti¬ 
mal gap and pressure [17,18]; between charge loss in the 
drift/conversion region and field focusing [4,19]; and be¬ 
tween high signal to noise ratio and low photon-feedback 
[20,7]. Ideally, considering the readout alone, one would 
wish to operate at the lowest possible drift field and high¬ 
est possible gain before the probability of feedback (or 
breakdown), recombination or attachment become impor¬ 
tant, as well as using a gas mixture and geometry that 
show a minimal sensitivity to the geometrical tolerances, 
[17]. For microbulk Micromegas, these aspects have not 
been discussed in full yet. 

We focus in this work on the modelling of microbulk 
Micromegas operated under Xenon and trimethylamine 
(TMA) for various admixtures and pressures. We benefited 
from the detailed experimental survey performed in [4], 
and the microscopic code Garfield-|--|- [21] in conjunction 
with the field solver Elmer [22] and Gmsh [23]. Considering 
that the Penning and recombination characteristics of the 
mixture might be at the core of an intriguing idea for direc¬ 
tional Dark Matter detection [24,25] and provided the only 
existing estimate for those properties is approximate [26], 
we attempted to extend the results from the microscopic 
electron tracking to include those additional effects. 

2. Simulation 

We followed for the simulation of the Micromegas 
response an approach based on the newly developed 
Garfield-|--|- [21]. The framework allows to track electrons 
by considering their elementary interactions with atoms 
and molecules, through the electron cross-sections tabu¬ 
lated in Magboltz [27]. It provides a more generic way of 
evaluating electron trajectories than the implicit thermal- 
ization ansatz of the hydrodinamic formalism (also avail¬ 
able in Garfield) while at the same time it allows accessing 
the fluctuations in the avalanche multiplication process. 
The latter are essential for a proper description of the 
energy resolution in gaseous detectors, that is attempted 
here. Ability for describing the avalanche fluctuations [5,6], 
and electron transparency through micro-meshes [28] has 
been demonstrated with this approach earlier. Nonetheless, 
the simulation framework does not contain at the moment 
libraries for the computation of electron-ion recombina¬ 


tion or atom-atom/molecule-atom transfer reactions, that 
are of relevance for the present work, and that have been 
hence evaluated externally. 

For the 3D-field calculation we relied on the open 
source packages Elmer and Gmsh, that are interfaced with 
Garfield-|--|-. The mesh was generated through tetrahe¬ 
drons on a elementary cell, that was then replicated to 
reproduce a volume large enough to avoid fringe effects. 

3. Micromegas geometry and electron transmission 

3.1. Geometry 

The geometry of the Micromegas chosen for this work 
is the one used in [4], an upper view of which is shown 
in Fig. 1, taken with a microscope. Its pattern consists of 
holes of diameter ^ = 48 ^m ± 2 /xm(sys) ± 2 /rm(sta) 
arranged as the vertexes of equilateral triangles of sides 
p = 100 pm ± 2 /im(sta). The systematic uncertainty in 
the hole diameter comes from its slight eccentricity, so the 
given range actually comprises the smallest and largest cir¬ 
cles compatible with the hole boundaries (dashed circles in 
Fig. 1). The statistical uncertainty comes from the measur¬ 
ing procedure. We have chosen in the following (jj = 50 pm 
and p = 100 pm, that were in fact the target values pro¬ 
vided to the photo-lithographic mask. The 3D electric field 
modelling was performed in the elementary cell shown in 
Fig. 2, that extends along the drift/conversion region over 
1 cm height. The thickness of the copper cladding (nano- 
coated with gold) was 5 ^m and the multiplication channel 
along the polyimide core was emptied in simulation by 20 
additional pm. The actual assumption is not critical for the 
simulated hole gain (representing a 20% effect at most) and 
is immaterial for the computation of the electron transmis¬ 
sion probability. The idea here is to eliminate any charge 
being lost to the insulator through diffusion, an effect that 
would likely cause gain transients and field reconfiguration 
(see [29] and references therein), and that is known to be 
small in microbulk Micromegas. The actual amount of poly¬ 
imide that is chemically etched inwards from the channel 
walls is not precisely known. By replicating the field in the 
elementary cell, the simulation was extended to a region of 
approximately 10x10 holes, ensuring no charge loss outside 
the geometrical boundaries under the conditions studied. 

3.2. Relative fraction of collected charge 

During Micromegas operation, it is difficult to unam¬ 
biguously isolate the fraction of charge directly gained in 
amplification from the one lost in the drift region or in the 
drift-to-hole electron transmission. An experimental deter¬ 
mination of the optimal charge collection conditions for op¬ 
eration usually involves a gain scan as a function of the drift 
Held (and normalized to the maximum gain obtained along 
the scan), hereafter dubbed IF. In simulation, on the other 
hand, the ‘true’ electron transmission Fmc can be obtained 
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Fig. 1. Microscope image of the surface of the microbulk Micromegas 
whose performance has been simulated in this work. The surface 
characteristics are determined by the gold coating. The dashed circles 
represent the biggest and smallest circles compatible with the hole 
foot-print. The (red) continuous circle has the average diameter of 
both, yielding an estimate of ^ = 48/imdi2/im. 


out-of-hole diffusion statistically cancel out for such high 
electric field conditions. 

The simulated transmission and the relative charge frac¬ 
tion obtained experimentally are shown in Fig. 3 as a func¬ 
tion of the drift field (Ed), for gains around 200-1000 and 
for a TMA concentration of ~ 1%. A ^°®Cd source was 
used for the measurements, and the analysis performed by 
assesing the variations of the e = 22 keV peak. First and 
foremost, it is apparent that this type of comparison pro¬ 
vides an excellent cross-check tool for the assessment of the 
geometry and quality of the structure. ^ Electron transmis¬ 
sion is expected to depend on the geometry, electric field 
configuration and gas properties, the latter known to yield 
an adequate description of the drift velocity and diffusion 
coefficients in Xe-TMA [30]. Excellent agreement had been 
reported earlier, e.g., in [28], for the transmission through 
standard micro-meshes. 



Fig. 2. 3D mesh used for the field calculation (in color online). 

easily by determining the fraction of primary electrons en- 
tering the holes of the structure (i.e., crossing a plane de¬ 
fined by the upper copper layer in Eig. 2). 

Electrons were launched in simulation at 200 /rm from the 
holes’ plane, uniformly distributed on an elementary cell. 
The results presented were found to be independent from 
the injection plane and from the particular distribution 
chosen (as long as the injection distance was greater than 
200 fim). Pre-amplification in the drift region was found to 
be negligible even in the presence of sizeable Penning trans¬ 
fer rates, so during the determination of the transmission 
curves this process was neglected. The finite element mesh 
density chosen in Pig. 2 provided a good compromise in 
terms of computing time, with finer meshes not yielding any 
visible difference. The asymptotic behaviour of the electron 
transmission (before the onset of pre-amplification) agreed 
with the optical transparency of the Micromegas within 
5%, implying that field lines become largely perpendicular 
to the hole plane. Although realistic electron trajectories 
in gases involve an additional random motion, in-hole and 
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Fig. 3. Simulated electron transmission (up), and fraction of col¬ 
lected charge determined experimentally. The labels indicate the 
operating pressure in bar. From 1 to 10 bar the TMA concentra¬ 
tions are, respectively: 0.90%, 1.01%, 1.15%, 0.86%, 0.86%, 0.86%, 
0.79%, 0.79%, 0.79%, 0.79%. Simulation results are presented with 
super-imposed splines for guiding the eye. The horizontal line indi¬ 
cates the optical transparency. 


Despite the good agreement in Fig. 3, small deviations 
between measurements and simulation can still be found 
both at high and (specially) low fields. To the right of the 


^ Prior to this work, and due to a miscommunication, the hole pitch 
was thought to be 115/rm, resulting in a gross disagreement. 






















maximum of each series, the measured transmission is sys¬ 
tematically rightwards shifted, by about 10%. The effect 
exceeds any systematic error stemming from the calibra¬ 
tion of the power supply (3%), drift distance (2%), tem¬ 
perature during gas injection (< 3%), TMA concentration 
(10%) and TMA inelastic cross-sections (20%). Variations 
along these lines end short at explaining the observations, 
however when considering the assumed uncertainty in the 
reference geometry (2/im) this additional shift can be ac¬ 
commodated easily. 

To the left of the maximum of each series in Fig. 3, charge 
losses due to either attachment or diffusion outside the ac¬ 
tive volume could be a priori expected, translating into a 
dependence of the collected charge with the drift distance 
z. This seems to be incompatible, however, with the long 
electron life-times (^ 5ms) measured in the system [19,7], 
implying charge losses below 0.1% for the 1 cm drift em¬ 
ployed. Additionally, the X-ray collimation was performed 
down to a few mm^ at the center of a 3.5 cm diameter Mi- 
cromegas sensor, un unbridgeable distance for the emerg¬ 
ing photo-electron. Ballistic defect can be discarded down 
to a sub-% contribution, due to the 50/iS integration time 
and the known values for the diffusion, drift velocity and 
photo-electron cloud [7] . ^ The most obvious culprit of the 
decrease of the collected charge at low fields is the presence 
of losses due to charge recombination. This effect, studied 
systematically in [19], has been recently reported in an ion¬ 
ization chamber measurement under 30 keV X-rays [31], 
qualitatively agreeing with the behaviour reported here. 

Some ‘ad hoc’ arguments given earlier in [32] and based 
on the Onsager [33] and Jaffe [34] solutions, suggest that 
recombination in gases may be approximated by a steep 
function, characterized by: 

a, 

where Qoo is the charge collected in the absence of recom¬ 
bination, Qo /Qoo represents the fraction of charge that es¬ 
capes recombination at (near) zero field, and k is an ef¬ 
fective parameter that describes the steepness of the func¬ 
tion. This functional dependence can be fitted to the left 
part of each series in Fig. 3-down, with the help of an ad¬ 
ditional global constant that corrects for the arbitrary nor¬ 
malization introduced in experimental data. The fit param¬ 
eter 77.(0) = 1-Qo/Qoo is plot in Fig. 4 (circles), exhibiting 
the trend expected from recombination: high pressure and 
TMA concentration result in a reduction of diffusion (and 
in an additional increase of the ionization density in the 
first case) hence reducing the fraction of charge that escapes 
recombination at zero field. For 10 bar, e.g., this fraction 
does not reach 50%. Remarkably, as argued below, a low 
diffusion enhances the field focusing and electron transmis- 

^ Since the spatial dimensions of the primary ionization cloud are 
compressed at high pressure and the drift velocity is largely a function 
of E^/P, a pressure increase would reduce ballistic deficit and thus 
increase the collected charge for any given E^/P, in stark contrast 
to the observed behaviour. 


sion, thus the net effect on the collected charge during de¬ 
tector operation is much less severe, going from 3%(lbar) 
to 9%(10bar) at most (Fig. 4 squares). 
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Fig. 4. Fraction of charge being recombined for zero field, l-Qo/Qoo, 
(circles) and for the reduced drift fields [Ed/P) employed during 
detector operation (squares). 


Aiming at a more direct comparison, we performed a 
fit to the fraction of collected charge T through a 4- 
parameter model: 


T{Ed, Ea) 


(1 - n)T 

max[(l — 77)7] 


( 2 ) 


with T{Ed,Ea) = TMc{Ed/ 9 ,Ea), and g being a scaling 
factor that generally stayed within [0.9-1] except for 1 bar 
(0.81). The drift {Ed) and amplification {Ea) fields are de¬ 
fined as the voltage drop across the drift and amplification 
gap, respectively. 

The result of the fit is shown in Fig. 5, where the relative 
fraction of collected charge has been split in 2-halves. The 
left part shows the effect of recombination, that depends 
mainly on the reduced field {Ed/P) except for an up-down 
modulation following the increase-decrease of diffusion and 
the decrease-increase of the ionization density. On the right 
part, the transmission is plot as a function of the ratio of 
field strengths, as customary in experimental data. This 
representation is partly justified since in the absence of dif¬ 
fusion the electron transmission is largely determined by 
the density of electric field lines entering the holes, that is 
invariant upon a global scaling of both fields. In reality, and 
due to the presence of a gas, the right edge of the observed 
transmission plateaus varies within a factor x3, following 
the same trend as the transverse diffusion coefficient, D^. 
For a typical drift field Ed = 0.2 kV/cm/bar, simulated co¬ 
efficients range from dOO/im/vTm x Vh^ (1.5%TMA) to 
ISOO^m/yTm X Vbar (0.1%TMA), and these strong dif¬ 
ferences persist when the electron progresses further into 
the hole. Hence lower transverse diffusion leads to stronger 
field focusing, as intuitively expected. 

The consistent behaviour of E and the qualitative agree¬ 
ment with the analysis performed in [7], [19] and indepen- 
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Fig. 5. Relative fraction of collected charge T in data and for the 
model given in text (lines). Up: as a function of pressure, for concen¬ 
trations around 1%. Down: as a function of concentration, for 5bar. 

dently in [31] point to the presence of charge recombination 
in Xe-TMA at fields < 300 V/cm/bar, increasing with pres¬ 
sure and with the concentration of quencher. It is however 
impossible, within the available measurements and simula¬ 
tions, to assure that the observed recombination is entirely 
taking place in the conversion volume. A sizeable fraction 
of charge recombination could be taking place at the en¬ 
trance of the Micromegas holes. The increase in ionization 
density resulting from the increased pressure and quencher 
concentration can potentially have a similar impact there, 
leading to the recombination of the incoming electrons from 
the drift region with the outgoing ions produced in the pre¬ 
ceding avalanches. 


4. Gas amplification 

The microscopic tracking of the electron trajectories in¬ 
cluding avalanche multiplication is an intense task: when 
performed in a standard computer it leads to computing 
times up to lOmin x cpu per primary electron, for the high¬ 
est pressures and gains. Nonetheless, calculations for real¬ 
istic gains of several lOO’s as those needed for /3/30 exper¬ 
iments [7] become possible, and they were realized in the 
present work by occupying 10 cpu cores during approxi¬ 
mately 3 months. 

Normally, in order to describe the avalanche multiplica¬ 
tion in simulation it is necessary to determine the proba¬ 


bility (rp) of Penning transfer processes [38], that are dom¬ 
inated in the present case by the reaction: 

Xe* -b TMA ^ Xe -b TMA+ -b e" , (3) 

Such a probability was evaluated through the following 
procedure: i) we started from the (comparatively simpler) 
estimate performed earlier under a parallel plate approxi¬ 
mation [35], and obtained a reasonable description of the 
measured gain for any given amplification field after 3-4 
simulations; ii) an improved estimate was obtained from an 
interpolation of the simulated vs gain points to the mea¬ 
sured gain; iii) as customary, was assumed to be field- 
independent, therefore its final value and uncertainty was 
obtained from the statistical combination of the values ob¬ 
tained for the high and low field regions of each gain vs Ea 
series. In each case we compute the effective gain as: 

rh* = r(I - 7^) X m (4) 

The (< 10%) corrections due to transmission and recombi¬ 
nation were obtained from the fits in Fig. 5. Simulated and 
measured gain curves for the most representative cases are 
compiled in Fig. 6 for TMA concentrations around 1% at 
Ibar {vp = 0.25), 2bar (r^ = 0.23), 5bar (r^ = 0.12) and 
lObar {jp = 0.14). Figure 6-right shows the distribution of 
gains (m) around the average gain m, together with a fit 
to a Polya function in the form given in [6]. These distribu¬ 
tions play an important role on the achievable energy res¬ 
olution and maximum gain [5], and they are usually char¬ 
acterized through their width / = The distribution 

of avalanche gains becomes a nearly perfect exponential at 
high pressure (/ = 1), since operation at smaller reduced 
fields Ea/P is enforced [36,37]. This fact, together with the 
increased charge recombination, anticipates a deterioration 
of the energy resolution at high pressure as demonstrated 
below. 
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Fig. 6. Left: experimental systematics of the effective gain vs am¬ 
plification field with different labels indicating the operating pres¬ 
sure in bar. TMA concentrations are 1.22%(1), 1.01%(2), 1.24%(5), 
1.71%(10) and lines indicate the simulation results. Right: simulated 
gain spread and fit to Polya functions. 








4.1. Penning transfer probability 

The ability to capture the microscopic atom-molecule re¬ 
actions leading to the effective transfer probabilities build¬ 
ing up Cp has been demonstrated in [38] for single-wire and 
for parallel plate geometries, but only for Argon mixtures. 
It is hence reassuring to note the approximate agreement 
between the Penning trends obtained in this work and those 
coming from the modelling of the single-wire data obtained 
earlier by Ramsey and Agrawal in [39] (Fig. 7, stars). These 
latter values were extracted from the measured gain curves 
by following the same procedure employed in [38]. 



Fig. 7. Penning transfer probabilities obtained from the simulation of 
the multiplication process for different pressures and TMA concen¬ 
trations. Lines represent the simplified transfer model introduced in 
text, that assumes the dominance of the Xe* ^Pi state. Dashed lines 
provide a straightforward extrapolation of the transfer probabilities 
from the assumed 50pm amplification region to an infinite space. 


We attempted a simplified description of the overall Pen¬ 
ning systematics given in Fig. 7 by considering only the 
processes of excimer formation fceac, Penning transfer and 
quenching (fcp, kg), together with an effective lifetime of 
the state undergoing Penning, r*. The problem is further 
simplified by considering only the two lowest lying Xenon 
excited states Esp^ = 8.43 eV, Esp^ = 8.31 eV and neglect¬ 
ing Penning transfers from the Xenon excimers. The lat¬ 
ter possibility is a highly plausible one since their energy 
falls well below the vertical ionization potential of TMA: 
IP^ = 8.44 eV [40]. Under these assumptions, the Penning 
transfer probability can be simply written as: 


fp = fi 


h 


cPk, 


'P,l 


(1 - cyP'^k^^^i P cP{kp^i + kg^i) + 1 /t* 
_ cPkp^3 _ 

(1 - c)‘^P‘^kexc,3 + cP{kp^3 + kq^s) 


( 5 ) 


is labeled as /i, p. ^ Although a fit to eq. 5 requires only 
5 independent parameters, it was impossible to constrain 
kp ,3 from the available data, its value representing unfortu¬ 
nately an important contribution for the extrapolation to 
low concentrations. It is sensible, however, to consider only 
the energy transfers from the singlet state. Its energy is 
nowadays known to be compatible with the vertical ioniza¬ 
tion potential of TMA within 10 meV. This is a sub-thermal 
energy difference at ambient temperature and it represents 
a remarkable fine-tuning of Nature, supporting the specu¬ 
lation that a resonant transfer may be taking place. After 
this, only 3 independent parameters are left in the fit: 


^ (1 - cyP'^ai P cPa2 P as 

The result of the fit is shown in Fig. 7 (continuous lines). 


with parameters: 

ai = 0.0047 ±0.0004 bar-1 

02 = 3.1 ±0.1 (8) 

as = 0.010 ±0.004 bar (9) 

By making use of the known transfer rate of the singlet 
state fceac,i = 51 10® bar”^ s”! [41]: 

T* = =9.2±2.5ns (10) 

^3 ^ear,! 


Due to the strongly mo del-dependent analysis procedure, 
the approximate agreement with the life-time of the sin¬ 
glet state Ti = 4.3 ± 0.5ns [42] must be taken with care. 
Effective life-times of atomic states can increase easily by 
3 orders of magnitude for cm-scale gas cells due to photon¬ 
trapping [42,43], although the reduced nature of the ampli¬ 
fication hole (50/rm) seems compatible with the small en¬ 
hancement observed for rf in the present situation. In Fig. 
7 the result corresponding to t* —>■ oo has been overlaid, 
a situation that would correspond to the Penning trans¬ 
fer rates observable in an unconfined space (e.g. a typical 
drift/conversion region). 

4.2. Energy resolution for 22 keV X-rays 

Each primary X-ray releases about rie = s/Wj = 900 
primary electrons, given an average energy for creating an 
electron-ion pair of Wj = 24.8 eV (in pure Xenon). In the 
range of TMA concentrations studied here we consider Wi 
to be the same for both species, an assumption that should 
lead to an error on the estimate of the direct ionization of 
< 1 %. 

A priori, the energy resolution in present conditions can 
be expected to be dominated by the stochastic nature of the 
avalanche multiplication process / (Eig. 6-right), i.e. 91 ~ 
2.35\//lF//e. This is because, as suggested by [6], during 
the operation of a gaseous detector in practical conditions 


Sub-indexes 1, 3 refer to the singlet and triplet state, re¬ 
spectively, and the probability of populating each of them 


® These probabilities, if neglecting the atomic cascade, can be esti- 
mated in first order from Garfield and they approach /i ^3 ~ 0.5. 















this contribution always exceeds the intrinsic fluctuations 
of the initial number of electrons, determined by the Fano 
factor (e.g., F = 0.f5 ± 0.02 in pure Xenon [44]). In detail, 
a complete description requires yet additional fluctuations 
related to charge losses ((1 — TZ),T), signal to noise (S/N) 
and mechanical tolerances to be considered (see [7], for 
instance): 


^int = + / + TTT. + (1 — T) - 


^mech 


S/N 

He 


1 dm 

~ “JT 

m dtp 

ENC 1 

m He 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


where ENC is the equivalent noise charge at the amplifier 
input, and all other magnitudes have been defined. The 
impact that the variations of the hole size, characterized 
through Umech j have on the energy resolution (either those 
stemming from the deviations from an ideal circle or hole- 
to-hole size variations) is difficult to assess, so the proposed 
idealized treatment must be necessarily considered as ap¬ 
proximate. The estimate of Umech may be characterized by 
the magnitude and a typical diameter spread cr^, un¬ 
der the assumption that the cloud of 900 primary electrons 
is large enough to cover a statistically representative re¬ 
gion. Fig. 8 shows the gain dependence with the hole size. 
For small holes the field approaches more closely the par¬ 
allel plate limit, leading to an increase in gain. More im¬ 
portantly, the slope is nearly a factor x3 higher at high 
pressure (table 1 ). 


P [bar] 

gain 


1 

350 

1.1 ±0.3 

1 

1350 

1.3 ±0.4 

10 

100 

3.0 ±0.6 

10 

150 

3.7 ±0.8 


Table 1 

Relative gain variation with respect to the hole diameter for different 
pressures and gains. 


For operation close to the edge of the electron transmis¬ 
sion plateau it was observed in simulation that the concur¬ 
rent increase of transmission with the diameter size can lead 
to an approximately compensated situation for the effective 
gain 771*, representing a limited analogy with the classical 
compensation in standard Micromegas [17]. In this discus¬ 
sion we assume for simplicity that measurements were per¬ 
formed far enough from the transmission edge so that only 
amplification is modified by the hole diameter spread. 

Energy resolution measurements are shown in Fig. 9, to¬ 
gether with simulation results after including all identified 
contributions (continuous lines). The ENC has been ad¬ 
justed to describe the low-field (low-gain) region, providing 



Fig. 8. Gain as a function of the diameter size, relative to the ref¬ 
erence geometry {(j) = 50/im). Normalization obtained from the su¬ 
per-imposed linear fits. 

values in the range 2000-4000e“; a typical diameter spread 
of cr 0 = O.G/rm was assumed, with the remaining contribu¬ 
tions being taken directly from simulation and formula 11 . 
As shown, the contribution of the hole accuracy is strongly 
hinted by data under an assumption on the diameter spread 
that seems realistic. Indeed, the value for Umech indirectly 
derived for the large demonstrator in [7] is consistent with 
the present analysis, both for 1 bar and 10 bar data, if set¬ 
ting (70 = 1/im. The presence of topological domains with 
different average values of (j) (within the quoted 1 fim) due 
to the larger areas considered there ( 8 mmx 8 mm per pixel), 
could naturally explain the additional spread. 



Fig. 9. Energy resolution for different TMA concentrations and pres¬ 
sures 1.22%(lbar), 1.01%(2bar), 1.24%(5bar), 1.71%(10bar). Con¬ 
tinuous lines represent the results obtained with the full simulation 
described in text, while for the dashed lines the contribution coming 
from the variations of the holes’ diameter has been eliminated. 

The deterioration observed at high gains for P = 2bar 
points to the presence of feedback, an effect that can com¬ 
pletely dominate the detector response (most strikingly for 


















operation under pure Xenon, [20], where the observed res¬ 
olutions exceed by a factor x4 the ones obtained here). For 
the type of readout and gas mixture studied in this work, 
the energy resolution plateau for which sufficient S/N and 
low photon-feedback can be simultaneously achieved ex¬ 
tends along a comfortable field range of about 10-20%. 

The influence of each contribution to the energy resolu¬ 
tion can be more easily extracted from the effective vari¬ 
ance per electron, defined as: 

”■ = (215) "• 

and that is shown in Fig. 10 for three representative field 
scans. 

5. Discussion on Fano factor and beyond-intrinsic 
energy resolution in Xenon 



E [kV/cm] 

amp 


In the presence of Penning transfer reactions, the energy 
resolution is modified as: 


= 


-912 _ 2.352 


(F-Fp) 


n, 


e,p 


n, 


e,p 


ne,p .. Ngx / JP \ 
—^ = 1 + -^rp{Ed) 
He Ni 


(17) 

(18) 


It depends on the number of excitations (susceptible of Pen¬ 
ning transfer) relative to the number of ionizations that are 
generated by the primary particle, on the Fano factor 
of the Penning mixture, Fp, and on the transfer probabil¬ 
ity at the typical drift fields, rp{Ed). The determination of 
Fp is obscured in the present case by the hypothesized con¬ 
tribution Umechi that exceeds the naive upper value set by 
Fp = F for nearly all pressures. At 1 bar and 1%TMA ad¬ 
mixture, however, and in the range Ea = 45-50 kV/cm the 
Fano factor expectedly becomes the second contribution in 
importance just after the avalanche fluctuations (Fig. 10- 
up), while the Penning transfer becomes maximal (Fig. 7). 
In such conditions, a determination can be attempted. 

In order to extract Fp, a set of additional assumptions is 
needed: first we assume that the Penning transfer probabil¬ 
ity extracted from the modelling of the amplification pro¬ 
cess is similar to the one at typical drift fields, so rp{Ed) = 
rp{Ea) = Tp = 0.24 ± 0.04. Additionally, we assume that 
the fraction of excited states susceptible of Penning trans¬ 
fer is directly related to the number of scintillation pho¬ 
tons, under which with Wsc being the energy 

that it takes to create a photon in Xenon, Wsc = 60 ± 20 
eV [45-47]: 

^ = l-b-^rp = 1.10 ±0.04 (19) 

Ue Wsc 

If it is true (as speculated in the previous section) that the 
singlet Xe* state predominantly contributes to the transfer, 
becomes smaller and the ratio approaches unity. 
We will assume for simplicity that this situation is approx¬ 
imately comprised within the assigned uncertainty. A com- 
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Fig. 10. Variance per electron, v*, and the various contributions 
discussed in this work. Up: gain scan (Ibar). Middle: gain scan 
(lObar). Down: drift field scan (5bar). 
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parison with formulas 11 and 16 in the region Ea = 45- 
50kV/cm leads finally to: 


Fp = - f -rn-{l-r)- Omech - <JS/N 

Tie. 

= 0.20 ±0.06 


( 20 ) 


The evaluation has been performed for Ibar and 1% TMA 
and makes use of the experimental measurement of v* to¬ 
gether with the parameters estimated in the Micromegas 
model developed in previous sections. The quoted uncer¬ 
tainty is mainly of systematic origin, dominated by the es¬ 
timate of amech- 

The Fano factor determined through eq. 20 is compatible 
with the one for pure Xenon, and nearly compatible with 
the minimum expected value [48]: 

= ±(l-rp) =0.11 ±0.02 (21) 

The intrinsic energy resolution of the gas mixture is thus: 


91o,Xe-TMyl = 2.35 



(0.50%±0.08%> 


IMeV 

e 


( 22 ) 


to be compared with the one expected in pure Xe [44]: 


94o,Xe 



(0.45% ± 0.03%) 


IMeV 


(23) 


This intrinsic energy resolution of the mixture (labeled with 
subscript 0) is important in scenarios where the multiplica¬ 
tion process is not the limiting factor in the reconstructed 
energy (e.g. [15], [49]). From such a perspective, it is im¬ 
plied by this analysis that the energy resolution achievable 
in Xe-TMA admixtures will represent only a modest im¬ 
provement over pure Xenon, if at all. Since the assumed 
Penning transfer probability is already close to the asymp¬ 
totic value Tp ~ 0.30 (eq. 6), a strong improvement seems 
unlikely for virtually any Xe-TMA gas admixture. 

The analysis presented suggests that the quality of 
microbulk Micromegas and microscopic modelling has 
reached a situation where the Fano factor of a gas mixture 
and its intrinsic energy resolution may be determined from 
the energy resolution obtained with X-rays (at least for 
certain operating conditions). A more systematic exper¬ 
imental verification of the indirect procedure here intro¬ 
duced is still needed, however, normally involving the use 
of different wafers or positions within the same wafer. On 
the other hand, the steady growth of the electron counting 
techniques performed on similar amplification structures 
can provide the necessary complementary path [49]. 


6. Conclusions 


interpreted as the result of electron-ion recombination at 
low drift fields, and the loss of field-focusing at high drift 
Helds (modulated by the transverse diffusion). A proper 
description of the gas amplification requires the presence of 
Penning transfers, expected due to the proximity between 
the ionization potential of trimethylamine and the energy 
levels of the excited states in Xenon. The probability of 
such processes, reaching maximum values in the range 
Tp = 20-30%, can be interpreted with a simple model that 
involves the resonant energy transfer from Xenon atoms in 
the ^Pi state to trimethylamine molecules. 

The energy resolution shows a degradation at high pres¬ 
sure slightly (but systematically) beyond the expected con¬ 
tributions, and it has been shown to be compatible with a 
higher sensitivity of the amplification process to the me¬ 
chanical tolerances in such conditions. In that respect, the 
accuracy of the hole manufacturing process has been eval¬ 
uated to be in the range < 0.6 ^m (in a 1 mm^ region) 
to 1 /im (in a 100 mm^ region), values beyond which the 
simulation becomes incompatible with the measurements. 

At low pressure, where the Penning transfers are higher 
and the influence of the mechanical accuracy was estimated 
to be negligible, the present work yields a Fano factor in 
Xe/TMAofFp(lbar, 1%TMA) = 0.20±0.06, together with 
an increase of the primary ionization by 10% ± 4% relative 
to pure Xenon. The intrinsic resolution of the Xe-TMA gas 
mixture amounts to: 

94o,Xe-TMA = (0.50% ± 0.08%) X i/lMeV/£. 
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We have shown that the behaviour of Micromegas 
manufactured in the microbulk and operated under 
Xenon/trimethylamine mixtures is well described by state 
of the art microscopic simulations in a broad range of con¬ 
centrations and pressures. When irradiated with 22 keV 
X-rays, the fraction of charge arriving at the holes can be 
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